Net-baryon-, net-proton-, and net-charged particle kurtosis in heavy ion collisions 

within a relativistic transport approach 
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We explore the potential of net-baryon, net-proton and net-charge kurtosis measurements to inves- 
tigate the properties of hot and dense matter created in relativistic heavy ion collisions. Contrary to 
calculations in a grand canonical ensemble we explicitly take into account exact electric and baryon 
charge conservation on an event-by-event basis. This drastically limits the width of baryon fluctua- 
tions. A simple model to account for this is to assume Poisson distributions with a sharp cut-off at 
the tails. We present baseline predictions of the energy dependence of the net-baryon, net-proton 
and net-charge kurtosis for central (b < 2.75 fm) Pb+Pb/Au+Au collisions from Ei a b — 2A GeV 
to y^ijvjv = 200 GeV from the UrQMD model. While the net-charge kurtosis is compatible with 
values around zero, the net-baryon number decreases to large negative values with decreasing beam 
energy. The net-proton kurtosis becomes slightly negative only for low y / 3jvjv". 
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A major aim of todays heavy ion programs at RHIC, 
CERN and GSI/FAIR is to explore the properties of ex- 
tremely hot and dense matter. Over the last decade 
event-by-event fluctuations of observables have been pro- 
posed as a novel tool to in heavy ion collisions provide 
information about the properties of the created mat- 
ter Proper observables based on these fluctuations 
should signal the onset of deconfmement [1, [H, [3] , reflect- 
ing the change in the degrees of freedom, or the criti- 
cal point, showing a non- monotonous behaviour with the 
beam energy [f| Q . Most fluctuation measures proposed 
are related to quadratic moments of the event-by-event 
observables. Among them are charged particle ratio fluc- 
tuations, charge transfer fluctuations, baryon number- 
strangeness correlations and charge-strangeness correla- 
tions. Higher moments and susceptibilities, however, be- 
come interesting as they are more sensitive to critical 
phenomena 0] and experimentally accessible. Recent 
lattice calculations by the RBC-Bielefeld lattice QCD 
Collaboration Q reveal that the ratios of higher suscep- 
tibilities in net-baryon, net-charge and net-strangeness 
number are around unity in the low temperature phase, 
increase towards a maximum at the deconfmement tem- 
perature and drop to values of to 0.5 in the deconfmed 
phase. These studies arc limited to zero baryo-chemical 
potential, conditions which can be met at RHIC or LHC. 

In this letter we investigate the potential of the net- 
baryon, net-proton and net-charge kurtosis under exper- 
imental conditions and constraints explicitly taking into 
account electric and baryon charge conservation. More- 
over, we provide baseline predictions for the energy de- 
pendence of the kurtosis of these quantities. Here the 
net-charge and the net-proton kurtosis are of special ex- 
perimental interest. Since the isospin susceptibility re- 



mains finite at the critical point critical fluctuations in 
net-baryon number are reflected in the fluctuations of 
the net-proton number Q. Critical phenomena leading 
to a non-monotonous behaviour on top of this baseline 
would be a signal for fluctuations due to a phase transi- 
tion or a critical point in upcoming energy scan exper- 
iments at RHIC and the future Facility for Antiproton 
and Ion Research (FAIR) at GSI. The numerical simu- 
lations are done within the microscopic transport model 
UrQMD in its latest version 2.3 0, E0, El • For previous 
event-by-event studies within the same model, the reader 
is refered to [3 [3, El. 

The suggested observables are the quadratic and quar- 
tic susceptibilities of baryons (£>), protons (P) and 
charged particles (Q): 
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where (SN 4 ) = ((N - N) 4 ) and {SN 2 ) = ((N - N) 2 ) 
are the fourth and second central moments of the re- 
spective distribution. The ratio of the quartic to the 
quadratic susceptibilities represents the kurtosis K(8N) 
times the variance of the distribution, a well known sta- 
tistical quantity: 
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Analysing the effective kurtosis instead of the kurtosis 
itself has the advantage that it eliminates the 1/N be- 
haviour as expected from the central limit theorem and 
removes the explicit dependence on centrality. The kur- 
tosis of a Gauss distribution is zero, -Kcauss = 0. In 
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FIG. 1: Effective kurtosis as from a cut Poisson distribution. 
The cut parameter evolves with the expectation value fi like 
in equation (J5J. For increasing fi the effective kurtosis first 
decreases until the the cut is predominantly limited by N to ^ . 



FIG. 2: Effective kurtosis of the net-charge, net-proton and 
net-baryon number distributions as a function of the width of 
the rapidity window as calculated from UrQMD for central 
Pb+Pb reactions at 158A GeV. 



comparison a distribution with K > (leptokurtic) has 
more weight in the peak and tail regions, while a dis- 
tribution with K < (platykurtic) has more weight in 
the flanks. Thus, platykurtic distributions are generally 
broader but have a lower peak and more narrow tails. 

Let us first explore what to expect in a transport model 
study of heavy ion collisions. While the lattice QCD cal- 
culations are performed in the grand canonical ensemble, 
transport approaches reflect the micro canonical nature 
of the underlying single scattering events even in the ther- 
modynamic limit. The choice of the thermodynamic en- 
semble, however, affects the fluctuations, e.g. discussed 
in [16] . In a grand canonical ensemble fluctuations of con- 
served charges are Poisson (or Gauss) distributed and the 
net charges are conserved only on average, while in heavy 
ion collisions net-baryon number is conserved exactly in 
each event. The fluctuations of net-baryon number in 
a subsystem are of the order of the root of the average 
net-baryon number in this subsystem. To a first approx- 
imation the exact conservation of net-baryon number is 
manifested in a cut-off at the tails of the net-baryon dis- 
tributions. What does this imply for the kurtosis analysis 
of heavy ion data in comparison to lattice QCD calcula- 
tions? To provide a simple example we introduce a cut 
in the Poisson distribution for the net-baryon number 



P„(N,C)=M(»,C)e-^ (4) 

on an interval [fx — C, fi + C] where /i is the expecta- 
tion value of the original Poisson distribution, C > 
indicates the position of the tail cut and A/"(/x, C) is a 
factor to normalize the modified distribution to unity. 
To model a (semi-)realistic case we assume that fluctu- 
ations in net-baryon number N are restricted even in a 



comparibly small rapidity window. In contrast to sta- 
tistical fluctuations, an increase of the rapidity window 
(and therefore of the average net-baryon number in the 
window) does not lead to stronger fluctuations. On the 
contrary, when the upper limit of the total net-baryon 
number N to t is approached the distribution changes to 
a (5-function (K$ = 0) due to baryon number conserva- 
tion. To capture this feature we thus assume a simple 
relation for the cut parameter 

c = ^( 1 -(£;) 2 )- (5) 

This means, fluctuations are statistical for a small subsys- 
tem and vanish in the limit of full phase space coverage. 

Since the Poisson distribution is a discrete distribution 
on integers, the lower integer number of /j, — C and ji + C 
is taken for explicit calculations. 

The resulting shape of the effective kurtosis is plotted 
in Fig. [T] with a — 3 and N to t = 416. It can be seen 
that for this choice the effective kurtosis of the cut Pois- 
son distribution is always negative and decreases with 
increasing fi until the cut affects the flanks of the distri- 
bution and weight is predominantly shifted to the peak 
leading to an increase of the effective kurtosis. The same 
qualitative picture can be drawn for a Gamma distribu- 
tion and the kurtosis of a cut Gauss distribution is always 
negative. The exact numerical values of the kurtosis cru- 
cially depend on the cut and finer details of the respective 
distribution. 

To study the quantitative behaviour let us now turn 
to the UrQMD investigation of the effective kurtosis. 
Here we vary the average net-baryon number directly by 
changing the size of the rapidity window at fixed ener- 
gies. In Fig. [21 we show the effective kurtosis of the 
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FIG. 3: Effective kurtosis for the net-charge, net-proton and 
net-baryon number distributions at midrapidity (\y\ < 0.5) as 
calculated from UrQMD at various beam energies for central 
Pb+Pb/ Au+Au reactions. The net-charge number at higher 
energies fluctuates severely due to limited statistics and is 
removed from this plot for clarity. 



FIG. 4: Effective kurtosis for net-charge, net-proton and net- 
baryon number distributions in acceptance windows around 
midrapidity covering an average net-baryon number of 150 
at each energy, as calculated from UrQMD at various beam 
energies for central Pb+Pb reactions. 



net-charge, net-baryon and net-proton number distribu- 
tions as a function of the rapidity window size in central 
Pb+Pb reactions at 158 A GeV. The latter is expressed by 
the mean net-baryon number in the acceptance, (N B a T B ). 
The acceptance windows are symmetric around midra- 
pidity. Broadening the rapidity window and thus going 
to larger mean net-baryon number leads to a decrease of 
the effective kurtosis until the size comes close to a 4-7T 
acceptance (in which the net-baryon number is strictly 

conserved, N B Q ^ B = 416), and the kurtosis shows the ex- 
pected steep increase. It is remarkable how well the shape 
reflects the assumption of a cut Poisson distribution. 

The effect of baryon number conservation is also seen 
in the energy dependence of the effective kurtosis. When 
going to lower energies, the constant rapidity cut im- 
plies that a larger fraction of the total (conserved) baryon 
charge is observed. Fig. [3] shows the effective kurtosis of 
the net-charge, net-proton and net-baryon number dis- 
tributions at midrapidity (\y\ < 0.5) at a range of beam 
energies from AGS to RHIC for central Au+Au/Pb+Pb 
reactions. Here one clearly observes the negative values 
for the net-baryon number fluctuations decreasing even 
further for lower energies. At RHIC energies the net- 
proton and the net-baryon number effective kurtosis is 
compatible with values of the order 1 and therefore in 
accordance with the latest data from STAR [l~7j . 

If the same fraction of net baryons is chosen (e.g. an 
average number of 150 net baryons is found in the ac- 
ceptance \y\ < 0.6 at 40,4 GeV and \y\ < 0.85 at 158A 
GeV, respectively), the resulting values for the kurtosis 
are comparable for all investigated energies (see Fig. [4j. 

Since protons contribute a large amount to the baryon 
number we expect a similar but less pronounced effect 



in their distributions. Only for lower AGS energies in 
Fig. [3] the effective net-proton kurtosis becomes slightly 
negative resembling the constraint on net-baryon num- 
ber. Also in Fig.[2]the effective net-proton kurtosis slightly 
follows the trend of the effective net-baryon kurtosis. Net 
charge is also a conserved quantity in a 4-7T acceptance but 
the fluctuations are less affected by this constraint. This 
is due to the large number of pions in a fixed rapidity win- 
dow that do not contribute to the average net charge but 
crucially determine the variance. Therefore net-charge 
fluctuations are less restricted by the constraint on the 
average net charge, which is determined by the protons. 
In our picture this is reflected in a less restrictive cut 
of the distribution. Even close to 47r charge fluctuations 
are still not completely supressed. Therefore, the effec- 
tive net-charge kurtosis is essentially not influenced by 
the conservation constraint, neither in Fig. [5] nor (within 
the large error bars) in Fig. [3J 

In RQMD QJ] and HIJING [H calculations at these 
energies we have also found negative values for the effec- 
tive kurtosis of the same order as are observed in UrQMD 
calculations. 

We have studied the energy dependence of the effective 
kurtosis of net-baryon, net-proton and net-charge num- 
ber within a transport model of heavy ion collisions that 
conserves charge and baryon number explicitely in each 
event. The net-baryon kurtosis was found to be negative 
for all energies. This observation, which is surprising in 
view of the lattice results given in [8j was explained by the 
exact baryon number conservation in heavy ion collisions. 
A variation of the rapidity window supported this inter- 
pretation and also a simple cut in the tails of a Poisson 
distribution can reflect the constraint on baryon number 
fluctuations qualitatively. A negative kurtosis shows that 
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moderate fluctuations are enhanced while the effect of 
large fluctuations is suppressed. A realistic treatment of 
heavy ion collisions at low ^/snn and, thus, large baryo- 
chemical potential yields drastic differences as compared 
to the lattice QCD calculations with respect to the net- 
baryon kurtosis. This is especially worrisome, because 
the onset of deconfinement is most probably located in 
this energy region. In our investigation, the effect of 
baryon number conservation alone is larger than the dif- 
ference stemming from a (potential) change of the un- 
derlying degrees of freedom as suggested by lattice QCD. 
This is, however, different for the net-proton kurtosis and 
its potential to be a signal of the critical point and the 
first order phase transition. Down to low SPS energies 
the background seen in the effective net-proton kurto- 
sis is only of the order ±1 and thus small compared to 
the expected signal from critical fluctuations [7j. How- 
ever, there are more possible effects in the experimental 
setup that might influence the kurtosis and support the 



trend to negative values, such as the selection of cen- 
trality classes and the mixing of distributions for single 
events. These influences will be subject to further stud- 
ies. The upcoming experiments of energy scan programs 
at RHIC and at FAIR will show if non-monotonous devi- 
ations from this baseline can be extracted signalling the 
phase transition or the critical point of QCD. 
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